Correlations in a band insulator 
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We study a model of a covalent band insulator with on-site Coulomb repulsion at half-filling using 
dynamical mean-field theory. Upon increasing the interaction strength the system undergoes a dis- 
continuous transition from a correlated band insulator to a Mott insulator with hysteretic behavior at 
low temperatures. Increasing the temperature in the band insulator close to the insulator-insulator 
transition we find a crossover to a Mott insulator at elevated temperatures. Remarkably, correlations 
decrease the energy gap in the correlated band insulator. The gap renormalization can be traced to 
the low-frequency behavior of the self-energy, analogously to the quasiparticle renormalization in a 
Fermi liquid. While the uncorrelated band insulator is characterized by a single gap for both charge 
and spin excitations, the spin gap is smaller than the charge gap in the correlated system. 



PACS numbers: 71.10.Fd, 71.10.Hf, 71.27. +a 
I. INTRODUCTION 

The role of electron-electron (e-e) interactions in solids 
is one of the central problems of condensed matter 
physics. The Hubbard model with local e-e interaction 
has become a paradigm for the description of electronic 
correlations in narrow-band materials. It has been used 
to investigate electronic correlations in metals and to 
study the correlation-driven metal- insulator transition. 1 
Much less attention has been paid to electronic correla- 
tions in band insulators (BI), since the lack of low-energy 
excitations rendered them less interesting. 

However, the discoveries of the quantum Hall effect 
and Kondo insulators showed that Bis are far from triv- 
ial, and recent progress in the topological classification 
of Bis 2 - demonstrates that our understanding of the in- 
sulating state is indeed incomplete. The quest for mate- 
rials with topologically non-trivial electronic structures 
suggests to explore heavier elements with strong spin- 
orbit coupling involving d or f electrons^ and raises the 
question about the role of electronic correlations. The 
common feature of these materials is that the constitut- 
ing atoms have partially filled shells and the gap - a 
hybridization gap - opens due to a particular pattern 
of inter-atomic hopping integrals. It has been proposed 
that similar characteristics apply to materials such as 
FeSi, FcSb2 or CoTiSb^ some of which exhibit strongly 
temperature dependent magnetic and transport proper- 
ties reminiscent of Kondo insulators. We call this class 
of Bis covalent insulators (CI) . 

Recently the evolution of a band insulator (BI) into a 
Mott insulator upon increasing the interaction strength 
has been studied in the context of the ionic Hub- 
bard model^&^ i 8 ' 9 ' 10 ! 11 : 12 a two-band Hubbard model 
with crystal-field splitting^ and a bilayer model with 
two identical Hubbard planes coupled by single-particle 
hopping j 1 16 i 17 i 18 Different scenarios have emerged 
with the possibility of an intervening phase and contin- 
uous or discontinuous transitions at critical interaction 
strengths. 



In order to study the properties of CIs with local e-e 
interaction we employ the dynamical mean-field theory 
(DMFT) ,!9i20i2ii22i23 In particu i ar we are interested in 

the nature of the interaction-driven transition from a co- 
valent to a Mott insulator, the possible existence of an 
intervening metallic phase, the evolution of charge and 
spin gaps and the single-particle self-energy as a func- 
tion of the interaction strength U. 

This paper is organized as follows: In Sec. [IT] we define 
the model and the methods chosen to study correlations 
in the covalent band insulator. The subsequent investi- 
gation is guided by the following questions: (i) How does 
the band insulator at weak coupling evolve into the Mott 
insulator at strong coupling (Sec. IIII Aj) ? (ii) What is 
the effect of correlations on the spectral function, and 
what happens when the temperature is increased in the 
correlated system (Sec. IIII Bj) ? (iii) Can we quantify cor- 
relation effects in a band insulator by means of concepts 
analogous to Fermi liquid theory (Sec. IIII C|) ? (iv) Re- 
garding the characterization of a simple band insulator 
as an insulator with identical charge and spin excitation 
gaps: Is this picture modified by correlations (Sec. IIII Dj) ? 
Our results are summarized and conclusions are drawn 
in Sec. UVl 



II. MODEL AND METHODS 

As a covalent insulator we denote a band insulator 
with partially filled identical local orbitals. This defi- 
nition implies that the band gap is a hybridization gap 
arising from a particular pattern of hopping integrals. 
Realizations of the covalent insulator include dimerized 
or bilayer lattices ; i 1 ' 16 ' 17 ' 18 quantum Hall systems with 
filled Landau levels or Haldane's model 2 ^ and the related 
model of Kane and Mele 2 - describing electrons on a hon- 
eycomb lattice with broken time-reversal invariance. We 
use a simple particle-hole symmetric model at half-filling 
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described by the Hamiltonian 

ka V fcCT y 

H(fe) = ( * _^ fe ) , (2) 

with two semi-circular electronic bands of widths 4t 
(t = 1 in the following) and dispersions efc and — efc, re- 
spectively, corresponding to two sublattices coupled by 
the fc-independent hybridization V and a local e-e in- 
teraction of strength U. Here n,i aa = a\ a a i(y measures 
the number of electrons with spin a =|, j on site i of 
sub-lattice a — a, b. 

We use the DMFT approximation to calculate the lo- 
cal single-particle propagator and the local spin suscep- 
tibility, quantities which within DMFT depend on the 
lattice only through the non-interacting density of states 
and thus are independent of a particular realization of 
the CI. The single-particle self-energy obtained by 

DMFT is local and fulfils the equations 

G(iw n )I=^ ((iw„ + ,i - S(iw„)) I -H(fc))- 1 , 

k 

Gq 1 (iw„) = G- 1 Qun)-'SQu)n), (3) 

where I denotes a 2 x 2 unit matrix. The functional de- 
pendence of E(iw n )[Go, U], defined on the discrete set of 
Matsubara frequencies uj n = (2n + l)irT, on Go and U 
is determined by an auxiliary Anderson impurity prob- 
lem, for a solution of which we employ the continuous- 
time quantum Monte-Carlo (CT-QMC) algorithm^ For 
quantities obtained in the imaginary time domain, i.e. 
the spectral function and the dynamical susceptibility, 
the analytic continuation to the real-frequency axis is 
performed using the maximum-entropy method^ 



III. RESULTS AND DISCUSSION 

A. Phase diagram 

The non-interacting ground state (U =0) of our model 
is characterized by a gap in the spectral function of size 
A = 2V at any V > 0. In the V = limit we have two 
decoupled copies of a single-band Hubbard model with 
semi-circular density of states, a problem which has been 
extensively studied within DMFT . 23 ' 27 It is well known 
that upon increasing the interaction strength U at fi- 
nite, low temperature the paramagnetic phase undergoes 
a discontinuous transition from a metal to a MI with a 
hysteresis in the interval U Cl < U < U C2 .— At high tem- 
peratures the hysteretic behavior is replaced by a contin- 
uous crossover. 
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FIG. 1: T-U phase diagram at fixed V = 0.5. The criti- 
cal values of the interaction strength U are determined from 
the double occupancy (see inset). Below a critical tempera- 
ture both band and Mott insulating solutions of the DMFT 
equations are found depending on the initial guess for the 
self-energy. For temperatures above the critical end point of 
the coexistence region, there is a regime where the spectral 
function has a single peak at the Fermi energy accompanied 
by broad Hubbard bands (see Fig. [3]). Inset: The double oc- 
cupancy D — (ji^ni) as a function of U for T = 1/30 (circles) 
and T — 1/60 (squares) indicates the phase transition from 
the correlated BI (larger D) to the MI (smaller D). Both 
values of the double occupancy are shown in the region where 
two solutions are found in the DMFT self-consistency cycle. 
Note that all energies are given in units of the hopping integral 
t = 1. 

In Fig. [1] we show the phase diagram in the T-U plane 
obtained for V — 0.5. At low temperatures two dis- 
tinct phases exist separated by a discontinuous transi- 
tion. The phase boundaries were obtained by calculating 
the double occupancy D = (nfnj,), shown in the inset 
for two selected temperatures. In a finite range of U 
values we find two stable self-consistent solutions of the 
DMFT equations.— The two phases, adiabatically con- 
nected to the band insulator (U = 0) and the Mott in- 
sulator (U — » oo), respectively, both exhibit a gap in 
the single-particle spectral function as discussed below. 
No signature of an intermittent metallic phase is found. 
The calculated phase diagram resembles that of a single 
band Hubbard model including the existence of a critical 
end point of the discontinuous phase transition, a weak 
T-dcpcndence of U c i and a considerable increase of U C 2 
upon lowering the temperature.— 



B. Single-particle spectral function 

The evolutions of the spectral function A(u>) = 
— Im G(uj + i0 + )/7r along a horizontal T = 1/30 and 
a vertical U = 5 line in the phase diagram of Fig. [T] are 
shown in Figs. [5] and [31 respectively. Remarkably, start- 



FIG. 2: Local spectral function A(lj) at fixed interlayer cou- 
pling V = 0.5 and temperature T = 1/30 for various values 
of the interaction strength U. The non-interacting density of 
states (U = 0) has a charge gap A c = 2V = 1. The gap in 
the correlated band insulator shrinks with increasing U until 
a discontinuous transition to a Mott insulator occurs with a 
hysteresis region 5.35 < U < 5.82. For U — 5.79 both the 
band (solid line) and Mott insulating (dashed line) solutions 
are displayed. All energies are given in units of the hopping 
integral t — 1. 

ing from U — the gap in the spectral function shrinks 
with increasing interaction strength U. At the same time 
the incoherent Hubbard bands evolve and spectral weight 
is transferred to them. The gap is well distinguishable 
throughout the entire interaction range except for a small 
region in the vicinity of U C2 where thermal broadening 
smears the strongly renormalizcd spectral features. 

At low T the spectral function in the band insulator 
phase close to the transition region consists of well distin- 
guishable low-energy quasiparticlc bands, separated by 
the hybridization gap, and incoherent Hubbard bands 
similar to the single-band Hubbard model. Increasing the 
temperature the spectral gap is filled in while the quasi- 
particle bands lose their spectral weight. At T = 1/10 the 
dip at chemical potential has vanished completely and 
a single peak remains. This peak smoothly disappears 
upon further increasing the temperature and at T = 1/5 
only two broad Hubbard bands remain in the spectrum 
reminiscent of the Mott insulator above the critical tem- 
perature of the metal-insulator transition^ In Ref. it 
was shown that this insulator-to-bad-metal crossover is 
reflected also in the dc and ac conductivity and is accom- 
panied by a substantial increase of the spin susceptibility 
which follows the Curie- Weiss law at high T. 



C. Gap renormalization and self-energy 

The spectral densities of Fig. [2] reveal a reduction of 
the charge gap in the CI phase with increasing U. In 
this section we analyze this behavior which is quantified 
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FIG. 3: Local spectral function A{ui) at fixed hybridization 
V — 0.5 and interaction strength U = 5 for various values of 
the temperature T. Upon increasing T the correlated band in- 
sulator (split peak around u — plus Hubbard bands) shows 
a crossover to a Mott insulator (broad Hubbard bands and 
reduced spectral weight near uj = at T = 1/5). At inter- 
mediate temperatures (T = 1/8 and T = 1/10) there is a 
single peak at the Fermi energy accompanied by broad Hub- 
bard bands. Note that the curves are shifted by a vertical 
offset for clarity. 

in Fig. [7] showing the charge gap A C (U) deduced from 
the spectral densities. In the following we derive the gap 
renormalization in two different ways from second order 
perturbation theory. Our derivation closely follows the 
approach of Ref. [29 . 



1. Renormalization of the charge gap deduced from the 
total energy 

The charge gap A c in the state with N particles is 
defined as 

A c = (E{N + 1) - E{N)) + (E(N - 1) - E(N)), (4) 

where E(N) is the ground-state energy of the system 
with N particles. The first correction to ([4]) from a per- 
turbative expansion in U comes from the 2nd order di- 
agram. Using the factorization of the joint density of 
states in the limit of infinite dimensions the second order 
correction to the ground-state energy can be written as 



Hi 

L 3 



E 

Pl,P2 ,P3,P4 



(1 - n PlT )n P2T (l - n P3l )n Pil 



c pi 



C P2 



1 I", 



(5) 

where n Pi is the occupation number and e Pi is the en- 
ergy of the single-particle state with index p.;. The non- 
interacting Af-particle ground state is a band insulator 
with all states with energies e < — V filled. The non- 
interacting (N + l)-particle ground state is obtained by 
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filling the lowest energy state of the empty conduction 
band e po = V (wc choose spin J. from the two possibili- 
ties). Keeping only the terms that do not vanish in the 
thermodynamic limit and using particle-hole symmetry, 
which requires that A c = 2(E(N+1)-E(N)), we obtain 
the second order correction to the charge gap 
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(0) 



Note that while £ (2) is an extensive quantity, the differ- 
ence in (|6|) remains finite when L — > oo. Introducing the 
single-particle density of states D and its Laplace trans- 
form F 



D(e) 



dee- Ae £)(e), (7) 



([B| can be rewritten ai 



,2!) 



A ( 2 ) = -2U 2 J dei de 2 de 3 J D(e 1 ) J D(e 2 )^(e 3 
x (1 - n 1 )n 2 



l-n 3 



€i - £ 2 + £3 - V ex ~ e 2 + V - e 3 



-4C7- 



poo 

/ dA sinh(AV)F 3 (A), 
Jo 



(8) 



where the fixed spin index has been dropped for simplic- 
ity. Introducing the bare gap A° = 2V we can write the 
renormalized charge gap as 



A, = A 



AW 



dA sinh 



A " 
A^ ) F 3 (A). 



(9) 



This expression shows that the reduction of the charge 
gap does not depend on the details of D(e), but rather 
on its overall characteristics such as the total bandwidth. 
If A° is small compared to the total bandwidth we can 
linearize the expression ^ to obtain 



A, = A? 1 - 2U 



dA F 3 (X)X 



(10) 



In this limit the gap acquires a simple multiplicative 
renormalization, which is closely related to the quasipar- 
ticle mass renormalization as shown below. The pertur- 
bative results are compared with the numerical data in 

Fig.m 

We close this section with two remarks. First, equa- 
tion ([5]) relies on the equality of the local and the total 
density of states, which is where the concept of a covalent 
insulator enters the algebra. Second, the physical origin 
of the reduction of the gap is best seen in ©. Adding 
a single electron to the insulating states blocks scatter- 
ing processes with a contribution of the order — 1/4V but 




FIG. 4: The self-energy within second-order weak-coupling 
perturbation theory at zero temperature (see Eq. (JTTJ) for 
V — 0.5. The real part (solid line) is linear around ui = 0; 
the slope (indicated by the dot-dashed line) determines the 
Z factor. The imaginary part (dashed line) is gapped for 
|w| < 3V. The dotted vertical lines indicate the bare gap 
A° = 1. The scale of the vertical axis is U 2 /t. 

adds the same number of processes contributing —1/2V 
and therefore leads to an overall gain in the correlation 
energy in the (N + l)-particle state and thus a reduction 
of the charge gap from its non-interacting value. 



2. Renormalization of the charge gap deduced from the 
self-energy 

An alternative derivation of the gap renormalization 
can be obtained from the pcrturbative calculation of the 
self-energy in the insulating ground state. Using the fac- 
torization of the joint density of states as in the previous 
section, the second order contribution to the self-energy 
is written as 

= U 2 Jdex de 2 de 3 D(e 1 )D(e 2 )D(e 3 ) 



(1 - ni)n 2 n 3 + ni(l - n 2 )(l - n 3 ) 



(11) 



(1 - £2 + U - £3 



where the spin indices were dropped as in the previous 
section. For — 3V < uj < 3V the denominator of the 
integrand remains negative throughout the entire inte- 
gration range, the self-energy is therefore real and can 
be expressed using the Laplace transform (J7J) of the non- 
interacting density of states 



EH 



-2U Z 



dA sinh(Aw)F 3 (A). 



[12) 



The renormalization of the band gap is obtained by 
searching for the pole f2 of the renormalized propagator 
of the lowest unoccupied state e po = V: 



n-v-s(n) = o. 



(13) 
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FIG. 5: Imaginary part of the self-energy as a function of 
Matsubara frequencies uj n for three values of the interaction 
strength U. In the correlated BI the self-energy has a negative 
slope at small ui n . In the coexistence region (U = 5.79) both 
solutions are displayed. 

In the small U limit we can replace E(f2) by S(F) and 
with A c = 2fl we recover equation ©. As in the previous 
section we can linearize (|12[) for small V in the interval 
M 5, V, which leads to 



where Z 



dRc S(w) 



duo 



O)=0 



(14) 



For | a; | > 31^ the second order self-energy acquires a fi- 
nite imaginary part and expression (JXSJ) is not applicable. 
The second order self-energy obtained by numerical inte- 
gration of (fTTj) is shown in Fig. [4] 

The linear behavior of Re £(u>) around the chemical 
potential is reminiscent of a Fermi liquid. It is not lim- 
ited to second order perturbation theory, but is a general 
consequence of a sufficiently fast vanishing of Im T,(u>) in 
the vicinity of the chemical potential and the Kramcrs- 
Kronig relations. While in Fermi liquids Im S(w) ~ u 1 , 
the existence of a gap in the spectral function of the co- 
valcnt insulator, as found in the numerical simulations, 
and the absence of a pole in £(u>) inside the gap, imply 
Im S(cj) = throughout the entire gap. As a result, 
Re of the CI phase closely resembles the self-energy 
of a Fermi liquid. This similarity is made evident in 
Fig. [5] where we show the self-energy on the discrete set of 
Matsubara frequencies: It is barely distinguishable from 
the analogous plot obtained for the single-band Hubbard 
model with comparable parameters^! 



D. Spin excitations in the band insulator 

Typically a band insulator is characterized by identical 
gaps for charge and spin excitations. So far we have iden- 




FIG. 6: Local dynamical susceptibilities in the correlated BI. 
At V = 0.5 and temperature T = 1/30 the imaginary parts 
of the spin susceptibility x. a (solid lines) and the bubble di- 
agram xb (dashed lines) are shown for U £ {2,3,3.5}. The 
susceptibilities are obtained from the QMC data by analytic 
continuation. The bubble diagram is calculated from the con- 
volution of the fully-dressed Green functions (see Eq. (|16|l )■ 
Thus the gap in Im Xs( CJ ) corresponds to the single-particle 
energy gap in A(uj), which is apparently larger than the spin 
gap observed in Im Xs( w ) m t ne correlated insulator. 

tified the insulator below the critical interaction strength 
as a "correlated band insulator" due to the fact that it is 
adiabatically connected to the band insulator at U = 0. 
Naturally the question arises whether spin and charge 
gaps remain indeed equal in the presence of correlations. 
To answer this question and to characterize more pre- 
cisely the correlated band insulator we evaluate the dy- 
namical spin susceptibility and compare spin and charge 
gaps. 

While the spectral function A(uj) is gapped in both 
band and Mott insulators (sec Fig. [5]), the spin gap 
is finite in the BI and zero in the MI. The spin ex- 
citation spectrum is reflected in the local dynamical 
spin susceptibility Xs(^)j which is calculated by a QMC 
measurement of the imaginary time correlation function 
Xs{t) = (S z (t)S z (0)) and the analytic continuation of 
its Matsubara transform to real frequencies. Here we use 
the maximum entropy method^ for the bosonic kernel 
according to 



1 f e~ TU} 
Xs{t) = - du i _ e _ fj J m x.(u). (15) 



In Fig.[S]the imaginary part of the spin susceptibility is 
shown on the real-frequency axis for V = 0.5, T = 1/30, 
and U S {2,3,3.5}. Similar to the charge gap, the spin 
gap also shrinks with increasing interaction strength. For 
a quantitative comparison of correlation effects on spin 
and charge excitations in the correlated band insulator 
Fig. [5] also shows the bubble diagram calculated from the 
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FIG. 7: Spin and charge gaps in the correlated BI as a func- 
tion of U for V = 0.5 as determined from the spectral function 
and the spin susceptibility at T = 1/30, respectively.— The 
thick solid line shows the charge gap obtained from second or- 
der perturbation theory (Eq. (0). The squares represent the 
Z factor (Eq. (|14|) ~). Inset: Discontinuous change of spin and 
charge gaps at the BI to MI transition. In the Mott insulator 
A s =0. 

convolution of the fully dressed local Green functions, 

In the non-interacting limit (not shown) Im xb(w) and 
Im Xs( w ) coincide and exhibit an energy gap of size 
2V = 1. In contrast at finite U the curves differ consider- 
ably from each other. In the spin susceptibility a promi- 
nent peak develops at energies lower than the charge gap. 
Furthermore spectral weight is suppressed at higher ener- 
gies which correspond to excitations from the split central 
peak to the Hubbard bands. The spin susceptibility is 
thus both qualitatively and quantitatively more strongly 
influenced by correlations in comparison to the bubble 
diagram, which contains correlation effects via the renor- 
malized propagator only. The effect of correlations on 
the band insulator therefore goes beyond the energy gap 
renormalization discussed in Sec. IIII CI 

The comparison of the spin susceptibility and the bub- 
ble diagram in Fig. [5] shows that charge and spin gap do 
not coincide in the correlated BI. For a quantitative anal- 
ysis we extract the gap values by a linear extrapolation to 
the frequency axis using the slope at the inflexion point 
of the spectral function and the imaginary part of the 
spin susceptibility for the charge and the spin gap, re- 
spectively. Fig. [7] shows the evolution of charge and spin 
gaps with increasing interaction strength in the corre- 
lated band insulator. As discussed in Sec. IIII CI the en- 
ergy gap renormalization at moderate coupling strengths 
is well described by the Z factor (Eq. ([13))) extracted 
from the slope of the self-energy at low frequencies. For 
U > 2.5 the spin gap is smaller than the charge gap, 
and A c — A s remains almost constant with increasing 



the interaction strength further. Therefore, in contrast 
to the non- interacting limit, the description of the band 
insulator as an insulator with identical energy gaps for 
both charge and spin excitations no longer holds once the 
interaction is strong enough. 

It is worthwhile to point out that different spin and 
charge gaps were also obtained in the half-filled one- 
dimensional ionic Hubbard model in the weakly corre- 
lated regime, when the ionic potential is smaller than 
the onsite interaction U As in the covalent insulator, 
A c and A s decrease with increasing U. However, the 
transition to the Mott insulator proceeds via an interme- 
diate insulating phase with bond order and a staggered 
modulation of the kinetic energy on neighboring bonds. 

Hints for a possible difference of spin and charge gaps 
as obtained in the above discussed correlated band insu- 
lator phase exist also from experiments on selected insu- 
lating materials. In their study of FcSi Schlcsinger et 
al.— pointed out the possible difference between spin 
and charge gaps in correlated insulators based on ev- 
idence for a larger charge gap in the Kondo insulator 
Ce3Bi4Pt3i22 While the single-particle charge gap can be 
measured with meV accuracy using photocmission spec- 
troscopy, the determination of the spin gap requires accu- 
rate susceptibility measurements to temperatures much 
lower than the gap energy. The relative gap difference 
for FeSi, e.g., is expected to be less than 30 % and could 
so far not be resolved from existing measurements.— 



IV. CONCLUSION 

We have studied correlation effects in a covalent band 
insulator using dynamical mean-field theory. In the ab- 
sence of correlations a band insulator is characterized by 
its band gap. A local Coulomb repulsion renormalizes the 
energy gap, which surprisingly shrinks when the interac- 
tion strength is increased. In second-order perturbation 
theory the gap shrinking can be traced to enhanced low- 
energy scattering phase space in the conduction band. By 
analogy to the quasiparticlc weight in Fermi-liquid theory 
a renormalization factor Z can also be introduced in in- 
teracting insulators based on the low-frequency behavior 
of the self-energy. The Z factor in the insulator deter- 
mines the energy gap renormalization. The simple one- 
gap picture of the band insulator breaks down for suffi- 
ciently large interaction strengths. In the correlated band 
insulator the spin gap is smaller than the charge gap. 
A discontinuous transition from the band to the Mott 
insulator occurs upon increasing the Coulomb repulsion 
at low but finite temperature. Close to the insulator- 
insulator transition the increase of temperature in the 
correlated band insulator with a split central peak and 
pronounced Hubbard bands leads to a crossover into a 
high-temperature Mott insulator phase with broad Hub- 
bard bands. 

The correlation-driven reduction of the energy gap 
is rare among the known materials. Nevertheless, this 
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mechanism provides a natural explanation of the uncom- 
mon gap overestimation in band-structure calculations of 
systems like FeSL— Tracing the difference between spin 
and charge gaps in the correlated band insulator to its 
physical origin remains a task for future work. 
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